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Abstract 

In the present work the zonal flow (ZF) growth rate in toroidal ion- 
temperature-gradient (ITG) mode turbulence including the effects of 
elongation is studied analytically. The scaling of the ZF growth with 
plasma parameters is examined for typical tokamak parameter values. 
The physical model used for the toroidal ITG driven mode is based on 
the ion continuity and ion temperature equations whereas the ZF evo- 
lution is described by the vorticity equation. The results indicate that 
a large ZF growth is found close to marginal stability and for peaked 
density profiles and these effects may be enhanced by elongation. 
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I Introduction 



The important role played by plasma flows for the reduction of energy trans- 
port in tokamak regimes of enhanced confinement is now widely acknowl- 
edged. Of particular interest for the improved confinement regimes are the 
nonlinearly self-generated zonal flows pQ. These are radially localized and 
strongly sheared flows in the poloidal direction. In experiments, zonal flow 
levels above the neoclassical prediction has been observed in connection with 
transport barriers [2] in the edge (H-mode barrier) and core (internal trans- 
port barriers) regions. Due to the strong shear stabilization of the driving 
instabilities associated with these flows, they are crucial for the dynamical 
self-regulation and saturation of the underlying turbulence and anomalous 
transport fluxes j3]. Theoretically, the generation and evolution of zonal flows 
due to drift mode turbulence has been extensively studied in recent years, 
both analytically [3] - [11] and in computer simulations using gyrokinetic 
[T2] - P3] and advanced fluid models [15J - [TT| . In particular, the gyroki- 
netic simulations of ITG mode turbulence reported in the Cyclone work Ref. 
[TS] indicate that there is a strong excitation of zonal flows close to marginal 
stability where the non-linearly generated flows were able to damp out the 
turbulence resulting in a non-linear up-shift in the critical temperature gra- 
dient needed to obtain transport for longer time scales. 

In the present paper, the excitation of zonal flows by toroidal ion-temperature 
gradient (ITG) turbulence is studied analytically. A system of equations is 
derived which describes the coupling between the background ITG turbu- 
lence and the zonal flow modes. The model used is an advanced fluid model 
for ITG modes [19] including effects of elongated flux surfaces. Previous 
work using this model indicates that the ITG growth is rather insensitive 
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to the shape of the flux surface [2TJ - [22]. The analytical technique used 
here closely follow Ref. [5] where the coupling between the zonal flow modes, 
driven by Reynolds stress forces, and the ITG turbulence is described by a 
kinetic equation for wave packets. The present work extends the previous 
study by studying effects of non circular geometry on the zonal flow growth 
rate. The result is an analytical expression for the growth rate of the zonal 
flow instability in elongated equilibria. 

The scaling of the zonal flow growth rate with plasma parameters is in- 
vestigated. In particular, the role of plasma elongation on the generation of 
zonal flow is examined. A resonance in the zonal flow excitation level is found 
close to marginal stability, consistent with the Cyclone work [18] . For peaked 
density profiles (small e n = 2L u /Lb), there is a substantially increased exci- 
tation of zonal flow when the elongation is increased whereas for most other 
cases the effects of elongation are weak. Moreover, the zonal flow excitation 
grows linearly with the wavenumber for zero collisional damping whereas for 
non-zero damping the zonal flow excitation is significantly reduced. 

The paper is organized as follows. In Section II the physical model for 
the toroidal ITG modes is presented. The equations describing the coupling 
between the background ITG turbulence and the zonal flow modes is pre- 
sented in Section III. Section IV is dedicated to the results and a discussion 
thereof. Finally there is a summary in section V. 

II Toroidal ion-temperature-gradient driven 
modes 

The description used for toroidal ITG driven modes consists of the ion con- 
tinuity and ion temperature equations. The effect of parallel ion motion is 
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weak in the reactive model and is hence neglected magnetic shear can, how- 
ever, modify the non-linear up shift as found in Ref. [10]. For simplicity, 
effects of electron trapping and finite beta effects are neglected in this work. 
The ion-temperature and ion-continuity equations can be written 

— + V • (riiVE + UiVu) + V • (riiV P i + ntf^) = (1) 
3 dT 1 

-n^ + n^V -Vi + V -qi = 0. (2) 

Here ve is the E x B velocity, is the diamagnetic drift velocity, v pi is the 
polarization drift velocity and v n i is the stress tensor drift velocity and qi, is 
the ion heat flux. The derivative is defined as d/dt = dj dt+p s c s zx V0- V and 
n, 0, Tj are the ion density, the electrostatic potential and the ion temper- 
ature, respectively. With the additional definitions n = 5n/no, <fi = e5(j)/T e , 
Ti = 5Ti/Tio as the normalized ion particle density, the electrostatic po- 
tential and the ion temperature, respectively. In the forthcoming equations 



r = Ti/T e , v± = p s c s y/L n , p s = c s /Vt ci where c s = ^T e /rrii, Vt ci = eB/rriiC. 
We also define Lj = — (dlnf / 'dr)~ l , rji = L n /L Tv e n = 2L n /R where R 
is the major radius and a, = r (1 The perturbed variables are nor- 

malized with the additional definitions n = L n / p s 5n/n , = L n J p s eS(p/T e , 
Ti = L n / p s STi/Tio as the normalized ion particle density, the electrostatic po- 
tential and the ion temperature, respectively. The perpendicular length scale 
and time are normalized to p s and L n /c s , respectively. The geometrical quan- 
tities are calculated in the strong ballooning limit (8 = 0, g (6 = 0, k) = l/n 
(Ref. [21]) where g{6) is defined by ojd{0) = uj^ngiO)). The equations 1 
and 2 can now be simplified to 



dn ( d 9 \ „ 2 7 d(j) d , ■ , . - 

-^7 - ^7 - oii— \7_l0 + e n £f— [(f) + t [n + Ti 

at \at ay I ay ay 



- [0, n] + [0, V10j + r [0, VI (n + %)\ (3) 
dfi 5 dfi ( 2\d4> 2 dfi 



dt d, Ttn9 dy + V l 3J dy 3dt 



- [0, Tj] + | [0, n] • (4) 

Here [A, B] = dA/dxdB/dy — dA/dydB/dx is the Poisson bracket. Lin- 
earizing the Equations 3 and 4 and using Boltzmann distributed electrons 
(hi =h e = 4>) gives 

(u - (1 - (1 + r) e n y) fcy + (u + (Xiky) kj^j + Te n gk y fi = (5) 
(u + ^re n #£;^ Ti - ( ~ 3) fc v + 3 w ) = ( 6 ) 

The corresponding dispersion relation and solutions are 

lOr 



= uj 2 



(l + k\) - uk y (l " (l + ~pj €n 9 - k l ( a i + 3^2) )(7) 

/ 7 5 5 \ 

+ T€ n gkl\n i -- + -(l + T)e n g + -a i kj_) (8) 

iTK)( 1 -( 1 + T)«- tl ("' + N) (9) 



k 1 

7 = x + y k2 V ^e n g fa - 77^) (10) 

where u; = cj r + ij and 

2 11 / 1 10\ 

^"3-27 + 4^ + ^U + 97j- 

FLR effects in the r] ith are neglected, however, they are important when the 
group velocities are to be determined. The group velocities (v g j = du r /dkj) 
are in the long wavelength limit (k\ « 1) given by, 

5r^ 



Vgx 



-k x k y (l + (1 + Vi) r - (l + y) e n g^j (12) 
Vgv = l( 1 -( 1 + ^-) e ng]. (13) 
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The instability of the mean flow in the linear regime sets in when the velocity 
of the mean flow modulations are close to the group velocities of the small 
scale wave packet. In Figure 1 the effects of elongation and e n on the r] r 
mode stability is illustrated. The growth rate (normalized the the electron 



diamagnetic drift frequency) as a function of rji with elongation («) and e n 
as parameters is displayed. The results are shown for r = 1 and e n = 0.2; 
K = 1.0 (diamonds) and k = 2.0 (squares) and e n = 2.0; k = 1.0 (+ curve) 
and k = 2.0 (* curve). 



Ill The wave kinetic equation and zonal flow 
evolution 

In this section the problem of how to construct the adiabatic invariant in 
ITG driven turbulence and deriving the evolution equations for zonal flows 
and ITG perturbations is summarized. The method has been described in 
detail in Ref. [5] (and References therein) and only a brief summary is given 
here. An alternative statistical approach, resulting in a modified wave kinetic 
equation (see Eq. 26 below), is presented in Ref. [TT] which also contains an 
extensive discussion of and comparison with the approach used here. In 
describing the large scale plasma flow dynamics it is assumed that there is 
a sufficient spectral gap between the small scale fluctuations and the large 
scale flow. The electrostatic potential is represented as a sum of fluctuating 
and mean quantities 

<f>{X,x,T,t) = $(X,T) + 4>(x,t) (14) 

where <£>(X, T) is the mean flow potential. The coordinates (X, T), (x,t) are 
the spatial and time coordinates for the mean flows and small scale fluctua- 
tions, respectively. From equations 3 and 4 we get after neglecting the FLR 
non-linearities, 

f + (1-(1 + r)e n9 ) f y -re n9 ^ = -[*,*] (15) 
^-3 re ^ + ^-3 (1 + r)e ^J^ = (16) 



Figure 1 
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Fig. 1 The ITG mode growth (normalized to the electron diamagnetic drift frequency) vs r\i 
with elongation and e n as parameters for r = 1 and e n = 0.2; n = 1.0 (diamonds) and n = 2.0 
(squares) and e n = 2.0; k = 1.0 (+ curve) and k = 2.0 (* curve). 
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Here, the interaction between the ITG perturbations have been omitted(see 
discussion after Eq. 26 for this.) In order to determine the generalized 

k = 

(pk + &kTi, where a k is to be calculated. Multiplying equation 16 by a k and 
adding it to equation 15 gives 



wave action density N k = \^k\ 2 we introduce the normal coordinates \I/ 



(J) k + a k f ik j + ^1 - (1 + t) e n g + a* (r)i - jj ^ 1 ' i L w j j q 
'1 \ dt. 



k 



- (^re n ga k + re n g^ = - $, <f) k + a k f ik (17) 
The normal coordinates are found if the equation is rewritten as in Ref. j5] 

where 

V k = l-(l+ T )e n g + a k (ri i -~{l + T)e n g\ (19) 

afc = — 3 r? (20) 

which gives 



-^(l-e n g + -^e n g)+%Je n g{r]i-r]ith) 

a k = " 2~7^ \ ( 21 ) 

Tfc - § (1 + t) e n g 

The linear relations between cf) k and enables one to express ^/ k and as 

2i^y k 

^k = 4>k + oi k T ik = - — — — (f) k (22) 



N » = l**r = TT7Sft 2 ( 23 ) 



i^ 2 



A fc = ^ f 1 - e n g + %e n g ) (24) 



2 V 3 

Ik = k y \fe n g (Vi -Vith) (25) 

The Eqs 22-23 and 25 describe the normal variables, the adiabatic invariant 
found from the normal variables and the linear ITG growth rate, respectively. 
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The wave kinetic equation see Refs. [5], [23] - [21] for the generalized wave 
action N k in the presence of mean plasma flow due to the interaction between 
mean flow small scale fluctuations is 

d M d , t _ s dN k (x,t) d (t _ s dN k (x,t) 

= 7 k N k (x,t) - AuN k (x,t) 2 (26) 

In this analysis it is assumed that the RHS is approximately zero (stationary 
turbulence). The role of non- linear interactions among the ITG fluctuations 
(here represented by a non-linear frequency shift Au) is to balance linear 
growth rate. In the case when 'y k N k (x,t) — AuN k (x,t) 2 = 0, the expansion 
of equation 26 is made under the assumption of small deviations from the 
equilibrium spectrum function; N k = N k + N k where N k evolves at the zonal 
flow time and space scale (Q, q x , q y — 0), as 

d 2 dN° 

-i(tt-q x v gx + ij k )N k = ky—$-^- (27) 

dk x it - q x v gx + ij k 
The evolution equations for the zonal flows is obtained after averaging the 
ion-continuity equation over the magnetic flux surface and over fast scales, 
including a damping term [25]. The average of the ion continuity equation 
over the magnetic surface and over fast small scales employing the quasi- 
neutrality the evolution of the mean flow is obtained 

|V>.« - rtft = (1 + r) VI (|*|*) + rVl (29, 

where it is assumed that only the small scale self interactions are the im- 
portant interactions in the RHS [26J. Using typical tokamak parameters 
(Ti = T e = Wkev, m = n e = 10 20 m- 3 , r = lm, R = 3m) y, = 0.78u H ^/{r/R) 
and Ua = 10~ 12 rii/T^ 2 and va is the ion-ion collision frequency, Tj is the ion 



temperature in electron volts. Using typical tokamak parameters it is found 
that n pd 50. Expressing the Reynolds stress terms in Eq. 29 in N k we obtain 

(-ifi - fxq 2 x ) $ = (1 + r + tS) J d 2 kk x k y \(p k \ 2 (30) 



where 5 is a k independent factor 

Vi-l(l + T)e n9 y (31) 



^ ^kky ( - 



Utilize equations 23, 28 and 30 gives, 

(-in - mt) = -nl (i + r + t8) 



d 2 kk 2 k x "^' k qxVgx (32) 
y x dk x Q - q x v gx + ijk 

where k x k y is substituted using the group velocity as 

it "gx 

X If / \ • 

1 + (1 + Vi ) r - (l + f ) e n 
It is now assumed that the short scale turbulence is close to marginal state 
(or stationary state, 7^ is small). Integrating by parts in k x and assuming a 
monochromatic wave packet N® = N 5 (k — k ) gives 

(n + ifxq 2 x ) (fi - q x v gx ) 2 = -q 2 x (l + r + t5) k 2 y N Q (33) 

In the special case of \x = 0, the third order dispersion relation for zonal flow 
Q reduces to 



^ — qx^qx "I - iqxk 



1 + T + r< y) ^-^N . (34) 



In expressing the zonal flow growth in dimensional form making use of the 
relation (A| + 7 jj!) / (4 7 |) iV = |0| 2 it is assumed that the mode coupling 
saturation level is reached [29J 

= -rV (35) 

CJ* KyLl n 
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IV Results and discussion 



An algebraic equation Eq. 33 describing the zonal flow growth rate including 
the effects of elongation is derived where the 3rd order dispersion relation is 
solved numerically and the zonal flow growth rates found are compared to 
the linear toroidal ITG growth rates. 

In Figure 2 the zonal flow growth rate (normalized to the ITG growth 
rate) as a function of e n (= 2L n /L B ) is displayed with elongation (k) as 
parameter for rji = 4, r = 1, q x — 0.3, fi = 0. The results are shown for k = 1 
(+ curve) and k = 2 (* curve). We note that a strong excitation of zonal 
flows (with 7/7/TG > 1) is obtained for small e n where a resonance is found. 
Similar results have been reported earlier in analytical calculations [10J and 
in the Cyclone non-linear (collisionless, electrostatic) numerical calculations 
of ITG mode turbulence [H]. In the numerical calculations there were a non- 
linear up-shift in the linear threshold from R/Lt ~ 4 to R/Lt ~ 6 and the 
other parameters were e n = 0.9, r = 1. The zonal flow growth is increased 
with increased elongation for small e n . 

Next, the effects of elongation and e n on the zonal flow growth rate is 
studied. Figure 3 shows the zonal flow growth rate (normalized to the ITG 
growth rate) as a function of k with e n as parameter and T]i = 4 and the other 
parameters as in Fig. 2. For peaked density profiles the zonal flow growth 
rate is increased with elongation and 7/7/TG > 1 is obtained whereas for flat 
density profiles (large e n ) the effects of elongation on the zonal flow growth 
are rather weak. 

Figure 4 illustrates the effect of collisional damping (//) and elongation on 
the zonal flow growth (normalized to the ITG growth rate). The normalized 
zonal flow growth rate is shown as a function of /x with parameter. 
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Figure 2 



Fig. 2 The zonal flow growth (normalized to the ITG mode growth) vs e n with elongation as 
parameter for r)i = 4, r = 1, q x p = 0.3, ji — and k = 1.0 (+ curve); n = 2.0 (* curve). 
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Figure 3 



Fig. 3 The zonal flow growth (normalized to the ITG mode growth) vs elongation (k) with e n 
as parameter for ^ = 4, t = 1, q x p = 0.3, ji = and e n = 0.2 (+ curve); e„ = 1.0 (* curve) and 
e n = 2.0 (squares). 
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The other parameters are as in Figure 3 with n — 1 (* curve) and n — 2 (+ 
curve). The zonal flow growth rate is decreased as the damping is increased 
and the effects of elongation are weak in this parameter regime. For /i — 50 
(typical value) a reduction of zonal flow growth with approximately 50 % is 
obtained. 

In Figure 5 the zonal flow growth rate as a function of the zonal flow 
wave number q x with elongation and damping as parameters is displayed. 
The other parameters are as in Fig. 3 with /j, — 0; k — 1 (* curve) and k — 2 
(diamonds) whereas /i — 50; k — 1 (+ curve) and k — 2 (squares). The effect 
of a non-zero damping is strong on the zonal flow growth rate and for zero 
damping the growth rate is linearly dependent on q x . Moreover, the effects 
of elongation are small. 

V Summary 

An algebraic equation which describes the zonal flow growth rate in the 
presence of toroidal ITG turbulence is derived and solved numerically in the 
strong ballooning limit. The model for the ITG driven mode is based on 
the ion-continuity and the ion-temperature equations. The present model is 
electrostatic and effects of trapping are neglected. The evolution of zonal 
flows is described by the vorticity equation and the time evolution of the 
adiabatic invariant in toroidal ITG turbulence is determined by the wave 
kinetic equation. This gives a system of equations that couples the zonal 
flow and the ITG driven mode perturbations. The equilibrium model used 
includes the effects of elongated flux surfaces. 

A strong generation of zonal flows is obtained for peaked density profiles. 
The effects of elongation on the generation of zonal flows for realistic tokamak 
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Figure 4 



2.5 



Fig. 4 The zonal flow growth (normalized to the ITG mode growth) vs fj, with elongation as 
parameter for rji — 4, r = 1, q x p = 0.3, e n = 1.0 and n = 1.0 (* curve); n = 2.0 (+ curve). 
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Figure 5 



Fig. 5 The zonal flow growth (normalized to the ITG mode growth) vs q x with elongation 
and p, as parameters for r]i = 4, t = 1, e n = 1.0 and p, — Ok — 1.0 (* curve) and k = 2.0 
(diamonds); \i = 50 re = 1.0 (+ curve) and re = 2.0 (squares). 
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parameters are rather weak, however, a resonance of the zonal flow generation 
is found close to marginal stability as in Ref. [TO] which is consistent with the 
results reported in Ref. [Hj. When damping is included in the model there is 
a significant reduction in the zonal flow generation for increasing zonal flow 
radial wave numbers whereas for zero damping the zonal flow generation is 
linearly dependent on the radial wavenumber. 

In short, this work indicates that for most parameter regimes, the effects 
of elongation on zonal flows driven by pure ITG modes are rather weak. 
However, the total effect on transport levels and turbulence may be significant 
since there is a resonance in the generation of zonal flows and this effect may 
be enhanced by elongation. For peaked density profiles, a strong excitation 
of zonal flows is found with 7/7/TG > 1- 
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